Magnetic field effect on the dielectric constant of glasses : an evidence of disorder 

within tunneling barriers ? 
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The magnetic field dependence of the low frequency dielectric constant e r (H) of a structural glass 
a — Si02+xC y H z was studied from 400 mK to 50 mK and for H up to 3 T. Measurement of both the 
real and the imaginary parts of e r is used to eliminate the difficult question of keeping constant the 
temperature of the sample while increasing H: a non-zero e r (H) dependence is reported in the same 
range as that one very recently reported on multicomponent glasses. In addition to the recently 
proposed explanation based on interactions, the reported e r (H) is interpreted quantitatively as a 
consequence of the disorder lying within the nanometric barriers of the elementary tunneling systems 
of the glass. 



PACS numbers: 61.43.Fs, 77.22.Gm, 72.20.Ht, 72.20.My 



I. INTRODUCTION 

Since the seventies pj, it is well established that at 
low temperature T, the properties of amorphous solids 
differ strongly from their crystaline counterparts. This 
is currently explained [2j within the tunneling two-level 
system (TLS) model, stating that the only relevant de- 
grees of freedom at low T are particles going back and 
forth between two neighboring sites of equilibrium. In- 
deed, the TLS model accounts for various properties of 
glasses below ~1K. A closer look, however, reveals that 
the TLS model often does not quantitatively account for 
the measured quantities. Since the late eighties, con- 
siderable efforts have been made to include interactions 
between TLS's in order to explain such discrepancies Q 
and justify the TLS model from first principles 0, [I| or, 
at least, to give a more consistent picture of TLS's 0, 
0- However, the lack of TLS microscopic information 
allowed to consider interactions either to be weak [j| or 
strong 0, 0j : these opposite starting points echo 
Leggett and Yu's question : "Is there anything really 
tunneling in glasses ?" 

The purpose of this work is to gain some microscopic 
information on TLS's by studying the magnetic field ef- 
fect on a glassy property, namely the 1 kHz dielectric 
constant e r . Indeed, for a charge Q, the tunneling trans- 
parency of a barrier of size a is strongly modified by mag- 
netic fields H of the order of H c — : this is due to 
the (H dependent) difference of quantum phase picked 
up by the various paths under the barrier Q- However, 
setting the usual TLS parameters, a ~ 0.2 nm and Q = e, 
with e the electronic charge, leads to H c ~ 10 4 T, a value 
so large that it seemed to explain Frossati et al's results 
[l0| reporting no detectable H effects on e r . Recently, 
however, detectable e r (H) effects were reported below 80 
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mK (and down to 1 mK) on a multicomponent glass y|, 
[TT| . Very strikingly, e r (H) is peaked around H c ~ 0.03 
T, a very low value accounted for by using interactions : 
at low T, it was argued that interactions correlate a large 
number N 3> 1 of TLS's and that the resulting spectrum 
resembles that of a single TLS with a charge Q = Ne 
(and a renormalizcd tunneling energy Ao). In this work, 
a detectable e r (H) effect is reported on a a — Si02+ x C y H z 
glass for 50 mK < T < 400 mK and H up to 3 T. Even if 
our glass is not really a model glass (such as the various 
forms of vitreous silica, e.g. Suprasil), the e r (H) depen- 
dence reported here shows that these effects are not re- 
stricted to the multicomponent glass studied previously 
0, favoring the idea that such effects are general 
in structural glasses. Furthermore our work shows, as 
stated previously |12| . that such measurements demand 
special experimental care since e r (H) effects are quite 
small: 5e' r (H)/e' r lies typically in the 10~ 6 — 10~ 5 range 
(e' r is the real part of e r ). This is why, after a brief de- 
scription of the experimental setup, it is explained how 
to overcome the difficult question of keeping T constant 
with great accuracy, while increasing H. The resulting 
t' r {H) effect is then accounted for by using numerical cal- 
culations of the H effect on the tunneling transparency 
of a disordered nanometric barrier. As explained below, 
this additional mechanism should not contradict the one 
based upon interactions. 



II. EXPERIMENTAL SETUP 

The sample was deposited on a a — S1O2 substrate 
as follows: i) a 0.1 /Ltm thick Cu layer was first evap- 
orated, ii) a L = 0.8 /xm thick a — Si02+ x C y H z layer 
was then deposited by a 13 MHz plasma vapor tech- 
nique using TetraEthylOrthoSilane mixed with He [li| 
iii) a 0.08 /nm thick Cu layer, followed by a 0.04 /im Au 
layer was evaporated. Both electrodes are 3 mm large 
10 mm long ribbons recovering each other along ~ 3 mm 
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(i.e., the capacitance surface is S ~ 9 mm 2 ) : the op- 
posite ends of the electrodes are free of glass deposit 
allowing an ohmic contact to be made with a Cu wire 
glued with Ag paste. These two Cu wires are soldered 
to the inner wires of cryogenic coaxial cables going from 
room temperature to the cold copper box embedding the 
sample : this box is related to the mixing chamber of 
the dilution fridge by a thermal impedance and contains 
a Ru02 resistive thermometer whose resistance will be 
called #r u o 2 throughout this work (to avoid any confu- 
sion with the parallel resistance R of the glass sample). 
The sample substrate was glued with varnish on a copper 
sample holder for thermalization. To decrease possible 
-ff-dependent heating in the electrodes, H was set par- 
allel to the electrodes (i.e., perpendicular to the 1 kHz 
measuring field £). Semi-rigid home-made coaxial ca- 
bles were installed at the top of the fridge and directly 
plugged in the Andeen-2500A capacitance bridge. The 
capacitance bridge was used in its parallel option, yield- 
ing C oc e' r and G oc e" but throughout this work it was 
chosen to report on C and R = 1/G due to the great 
similarity of the C(T) and R(T) curves (see Fig. 1) : 
this strongly confirms that the reported R at low T is 
not some parasitic edge effect but has the same origin 
that the reported C, namely the TLS's. 

From 300 K to 4.2 K, the capacitance C typically 
halves, while the parallel resistance R increases by a fac- 
tor ~ 100. Below 4.2 K, C decreases, reaches its min- 
imum for T rev , and then increases below T rev _ before 
reaching a ultra low T saturation value C sa t- Accord- 
ing to the standard TLS model the C decrease above 
T rev , is due to progressive freezing of the diagonal (or 
relaxational) part of the response, while the C increase 
below T rev . comes from the induced off-diagonal (or res- 
onant) part of the susceptibility : this effect enlarges as 
T decreases as do all quantum effects. 

III. OBTAINING Se' r (H) 
A. Proof that a non-zero Se' r (H) exists 

Figure 1 shows C(T) and R(T) curves when a 1 kHz 
field £ = W/0.8 = 3.75 MV/m is used. The high 
value of T rev _ ~ 0.3 K ~ 8T rev \£ — > 0) reveals that such 
fields correspond to the strongly nonlinear regime |14| . 
Such large fields were used with a large integration time 
(100 s) so as to reduce the bridge uncertainties c no i se 
(resp. r no i se ) over C (resp. R) measurements down to 
c-noise/C = ±10~ 7 (resp. r noise /R = ±10~ 5 ). Large 
fields were also used in [Tl| . 

The second key point to investigate small e r (H) effects 
is to be able to increase H while keeping T accurately 
constant. It demands either to set up a H- independent 
thermometry, which is extremely difficult over an ex- 
tended T range (it may have been the way used in @, 
fll|h or, at least, to correct T of the magnetoresistance of 
the thermometer. However, such a correction cannot be 




FIG. 1: C(T) (left axis) and R(T) (right axis) when a 1 kHz 
measuring field £ = 3.75 MV/m is applied to the glass sam- 
ple. Lower inset: 8r u o 2 (H) dependence measured assuming 
that the V(T) relationship does not depend on H (see text). 
Upper inset: Schematic effect of a negative RuC>2 magnetore- 
sistance : increasing H while keeping (9r u o 2 constant leads to 
a small decrease of T, i.e., ST(H) < 0. 



done with great accuracy since the standard technique for 
measuring the #r u o 2 (H) dependence is as follows [17| : i) 
For H = the relationship between the power V injected 
in the fridge and the resulting T is first recorded ; ii) For 
a given V, H is slowly raised to a given value, and the 
measured 6*r u o 2 {H) — 6*r u o 2 (0) at thermal equilibrium is 
attributed to the magnetoresistance of the thermometer 
: the corresponding results are shown in the lower inset 
of Fig. 1 and lead to corrections of the order of 2 — 5% 
on T when H is raised up to a few T. However, these 
0RuO 2 (H) measurements lie on the key assumption that 
the V(T) relationship does not depend at all on H : this 
cannot hold accurately, e.g., due to the dissipation of en- 
ergy resulting from the vibrations of the fridge within the 
magnetic field. These small and complicated effects can- 
not be ignored when seeking small e r (H) effects which, 
translated into a thermal equivalent, amount to an effec- 
tive variation of T lying in the 0.1 — 5% range (see T5)2) 
below). 

Going back to the C, R measurements, the fact that 
both V(T) and 0r u o 2 (? 1 ) relationships sligthly depend 
on H implies that T cannot be made accurately con- 
stant when H is increased. For practical reasons, the 
C(H),R(H) measurements were made letting V auto- 
matically adjust so as to keep the #r u o 2 value fixed at 
any H value : raising H while keeping #r u o 2 constant 
leads to a small shift 5T(H) of the temperature of the 
experiment as depicted in the upper inset of Fig. 1. 

The measured variations 6C(H) = \C{H) — 
C(0)] CO nstant # Ru o 2 are displayed in Fig. 2 for various 
T values measured at H = 0. The key point is that the 
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FIG. 2: SC(H) measured by keeping #r u o 2 constant while 
increasing H (the various T labeling the curves are at H — 0) . 
Inset: Corresponding data for SR(H) (same symbols). The 
fact that SR(H) and 6C(H) do not vanish at the same H 
cannot be explained by using only ST(H) effects but implies 
that a Se r (H) dependence exists in the sample. The error 



bars Cnoise = ±10 pF and r nois 
smaller than the data symbols. 



= ±10 b Gfi are much 



sions are accurate. Defining 5X(H) = X(H) — X(0) for 
any quantity X measured as a function of H keeping the 
value constant, one gets : 



'Ru0 2 



6T(H) + 6C mt .(H), 



(la) 



where the measured 6C(H) is shown on Fig. 2, 
SCi n t. (H) is proportionnal to the sought intrinsic vari- 
ation of Se' r (H) and throughout this work, all the partial 
derivatives with respect to T are taken at H = 0. Intro- 
ducing the corresponding quantity 5Ri n t.(H), an equa- 
tion similar to Eq. (la) is obtained for R measurements. 



SR(H) = (gf) ST(H) + SR mt . (H). (lb) 

Combining Eqs. (la) and (lb) allows to accurately 
eliminate ST(H), and to get : 



5C mt .(H)= v 




6R(H) - 6C(H) , (2a) 



with i] — ?— jrn an d a is defined by : 



8T(H) dependence cannot explain all the qualitative fea- 
tures of Fig. 2. For example, whatever ST(H) may be, if 
Se r (H) were zero, SC(H) and 5R(H) would vanish at the 
same H : this is clearly contradicted by the 50 mK, 100 
mK and 400 mK data reported on Fig. 2. This proves 
that a e r (H) dependence exists in our sample but that an 
accurate correction of the 8T(H) effects must be done on 
the data of Fig. 2. This correction cannot be done with 
the data of the inset of Fig. 1 since, as above evoked, the 
sought e r (H) is smaller than the accuracy of the 5T(H) 
estimate. 



B. Obtaining the order of magnitude of 5e' r (H) 

At this step, the existence of a detectable 5e' r (H) de- 
pendence is established but an accurate suppression of 
8T(H) effects is still missing. It is shown here that us- 
ing both 5R(H) and 8C(H) allows to draw from Fig. 2 
the order of magnitude of Se' r (H), i.e., to know its value 
within a scale factor 0.1 < n < 10 (see Eq. (2) below) : 
this interval might seem to be very large, but we will see 
in the next section that the physical interpretation only 
involves 111(77). 



1. Overview of the method 

The key idea is that since the relative variations of 
all the involved quantities are small, first order expan- 



SRjnt.(H) (fy) 

5C int .(H) (ff)' 



(2b) 



a is introduced here as an attempt to eliminate one of 
the three unknown quantities among ST(H), 5Ci n t.(H) 
and SRi n t.(H) since only two quantities, SC(H) and 
5R(H), are measured, a cannot be predicted since, for 
strong £ , a comprehensive treatment of the TLS dynam- 
ics, including quantum coherence effects, is still missing 
: at present, even the C(T) curve is not well accounted 
for |l4| . However, within the framework of the stan- 
dard TLS model, a and 77 cannot depend on H since, 
as we show below, Ao jma i is the only parameter of the 
TLS model which might slightly depend on H : neglect- 
ing, as above, any second order terms, both the real and 
the imaginary parts of Se r (H) (i.e., both SCi n t.(H) and 
5Rint.(H)) are expected to be proportionnal to the same 
small 5Ao )ma x(£0/Ao,maa! ( se e Eq. (4)). 

There is no similar argument stating that, generally, a 
and r) = l/(a — 1) cannot depend on T. It is shown in 
the section A of the Appendix, by using the shape of the 
curves of Fig. 2 and quite general arguments, that rj lies, 
at any reported T, in [0.1, 10], an interval which will turn 
out to be "small" enough to neglect any i](T) variations 
(see section IV). Basically, it is shown in the Appendix 
that the 5T(H) influence on the measured 5C(H) and 
5R(H) enlarges as H is increased : the quasilinear de- 
crease of 5C(H) 1 8R(H) seen in Fig. 2 at large H thus 
comes from ST(H) effects (at 400 mK it turns to a quasi- 
linear increase due to the fact that dC / dT > contrarily 
to the case of lower T data). 
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FIG. 3: SCint.(H) oc 5e' r drawn from Fig. 2 data by using 
Eq. (2a) with rj = 1. For clarity the error bar (see Eq. (3)) 
is reported only for one H at each T : it is much stronger at 
400 mK than below 100 mK where error bars hardly exceed 
the size of the data symbols. Inset: SC(H) measured just at 
T rev . (with a larger field £ = 8.75 MV/m to further increase 
the signal to noise ratio) : due to the dC/dT = at T rev . 
the use of Eq. (2a) is not necessary at low H (at higher 
H, ST(H) is no longer negligible and produces the SC(H) 
increase). Comparison of the data at T rev . to those at 400 
mK supports the idea that rj > at 400 mK. 



2. The most salient features of SCint. (H) 

Since 0.1 < rj < 10, the order of magnitude of 
SCi n t.(H) is obtained by setting rj — 1 and disregard- 
ing any possible rj(T) variations within this interval : the 
resulting transformation of the data of Fig. 2 through 
Eq. (2a) is shown on Fig. 3. 8Ci n t.{H) is positive for 
the lowest temperatures (except for three low fields at 
50 mK, see below) and negative at higher T, i.e., in the 
vicinity of T rev ,. 

The 200 mK data were not reported on Fig. 3 since 
they lie very close to zero. Indeed, it seemed more im- 
portant to report, in the inset of Fig. 3, on SC(H) mea- 
sured just at T rev , (measured with a larger field £ = 8.75 
MV/m to further increase the signal to noise ratio). For 
the inset of Fig. 3, at low H , the ST(H) contribution 
(see Eq. (la)) to 8C(H) is negligible since dC/dT = 
just at T rev , (at large H this is no longer true due to 
the increase of 5T(H), yielding the quasilinear increase 
of 8C(H)). The important point here is the decrease of 
5C(H, T rev ) at low H : it further confirms that SCi n t. (H) 
is negative in the vicinity of T rev , and that 77 > at 400 
mK. Note that the results of Fig. 3 were reproduced on 
a second similar sample. 

Let us translate the data of Fig. 3 to a thermal equiva- 
lent by using ST int .c = 5Ci nt .(H) / (remember that 
8Ci n t.{H) are calculated by setting 77 = 1) : it is found 
that 5T int .c(H) reaches at most 7% of T at 50 mK, 2% 



at 70 mK, 1% at 100 mK and 4% at 400 mK. Beyond 
chemical differencies between glasses, this might explain 
why, with an uncertainty on T of ±5%, previous studies 
reported undetectable e r (H) effects |l0j| . Adding quadrat- 
ically the uncertainties on both terms of the right hand- 
side of Eq. (2a), the final uncertainty c err about the 
SCint. (H) reported on Fig. 3 reads : 



V\ (Cnoise) 2 + {r r , 



8C_ dR, 

' dT dT' 



(3) 



On Fig. 3, c err was reported for each T : due to the 
strong decrease of dR/dT, it is much stronger at 0.4 K 
than at low T where it hardly exceeds the symbol size. 
However, even at 0.4 K, c err amounts to an uncertainty 
over T of ±0.5%, ten times smaller than in [lflj . 

Let us note that SCi n t. (H) / c err does not depend on the 
chosen r\ = 1 : the trends of 5Ci nt .(H) reported on Fig. 
3 are thus reliable at each T, even if the possible f]{T) 
dependence forbids any accurate comparison of SCi n t. (H) 
at different T . Comparing Fig. 3 with Strehlow's results 
[Tlj . our 8e' r (H) is somewhat smaller but lies in the same 
range. As in hj, our Se' r (H) might be peaked around 
H ~ 0.02 T for our lowest T, but this effect is too close 
to our c err to be systematically studied. 



IV. PHYSICAL INTERPRETATION 

Let us move to the interpretation of Fig. 3. First, note 
that, due to the T-dependence of the data of Fig. 2 and 
Fig. 3, the reported 5e' r (H) can be attributed neither 
to the coaxial cables nor to the "matrix" surrounding 
TLS's. Furthermore, the change of sign of 8Ci nt .{H) in 
the vicinity of T rev . excludes a role of the electrode-glass 
interfaces and further confirms that TLS's are responsible 
for the reported 8e' r (H). Last, it is shown in the section 
B of the Appendix that spin effects, if any, can be ruled 
out at low enough T, i.e., well below T rev , : this is why 
we will focus on the positive 5e' r (H) observed at T < 200 
mK. 



A. Going from 8C in t.{H) to 5A , ma x(H) 

For the part of the capacitance Ctls coming from 
TLS's, the usual result [14J below T rev , can thus be used 
: Ctls — C hi (E max /kBT) where ks is Boltzmann con- 
is the maximum inter- 



stant, E„ 
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max 



A 2 
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level spacing of TLS's, and the constant C ~ 0.03 pF 
is determined by fitting the increase of C(T) of Fig. 1 
below T rev , (details are given below). 

Since both A (see section B of the Appendix) and C 0] 
are if- independent, the reported Se' r (H) is interpreted as 
a Ao tmax (H) dependence, yielding : 
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15Ci n t\H) _ SA . max (H) 



C 



.max (0) 



ln( 



(H), 

Ao jm ax(0) 



(4) 



where the standard assumption A ma2; = Ao imaa ; is re- 
sponsible for the factor 2 in the left hand-side of Eq. (4). 

Let us note that the expression of Ctls used to derive 
Eq. (4) only holds accurately between 90 mK and 200 
mK since, at lower T, the C(T) curve tends to saturate. 
This failure of the fit at low T is due to the fact that the 
above expression of Ctls is well established only in the 
linear regime, and that a comprehensive treatment for the 
strong £ used here is still missing. However, this formula 
should capture the essential part of the physics since it 
expresses that, well below T rev _, the TLS susceptibility 
is due to the quantum transitions coherently induced by 
E. Since, physically, it is likely that tunneling plays a 
key role in the nonlinear Ctls at low T, Eq. (4) should 
finally hold, up to a multiplicative prefactor which can 
be disregarded just like r\ (see below). 



B. A possible microscopic explanation of 

A Ao Jmax (H) dependence means that the tunnel trans- 
parency is affected by H. This is well known in another 
kind of disordered insulators, namely Anderson insula- 
tors (e.g. lightly doped semiconductors) where electronic 
states are localized due to potential disorder. At very 
low frequency, transport in Anderson insulators occurs 
by hopping on a fractal percolating path strongly affected 
by H. H effects are smaller at high frequencies where the 
conducting path is broken since only the fastest hops con- 
tribute : note that, since they are disconnected from each 
other, each of these fast back-and-forth hops is a TLS. 
In both limits, H effects were accounted for ^]J, G3 by 
studying the tunnel transparency Ao between two sites 
separated by a disordered barrier depicted in the inset 
of Fig. 4 : tunneling between the two sites of the (qua- 
sisymmetric) TLS separated by the distance a is strongly 
affected by the potential disorder coming from the struc- 
tural disorder within the tunneling barrier itself. This is 
modeled by a three dimensional network of " impurities" 
(whose unit length is the elastic mean free path X e {) and 
Ao results from the coherent sum of all quantum paths 
along the impurity network : at H — 0, a path picks up 
a phase —1 each time it diffuses on an impurity whose 
energy e, is smaller than that of the TLS. Since the 
are drawn at random, Ao must be first averaged over 
disorder (noted hereafter with < >). 

Using this model to account for the measured Se' r (H) 
is obviously an oversimplification. Indeed, TLS's might 
be more complicated than a unique charge going back 
and forth between two sites separated by a: a TLS could 
be a group of charges, a 'molecule', and more generally 
any charged species with two spatially separated poten- 
tial minima. However, we show now that this model ac- 




1.E-02 



1.E-03. 



y— a = 800 A, e | 

a = 400A, P i 

^ — a = 100A? 



□ 100mKdataof Fig.3 



1.E-06 



1.E-05 |_| J |_| 1.E-04 



1.E-03 



FIG. 4: Inset: Schematic view of the disorder lying within 
the tunnel barrier of size ~ a between the two sites (gray cir- 
cles) of a TLS : the potential fluctuations are modeled by a set 
of "impurities" of interspacing A e ; whose energies are drawn 
at random either well above or well below that of the TLS. 
Main figure: Ao results from the coherent sum of all the paths 
linking the TLS sites : due to the additional quantum phase 
introduced by H for each path, < In Ao > increases with H, 
the H 2 dependence at very low H (unreported) turning con- 
tinuously to a H 5 behavior at higher fields (ensemble averag- 
ing over the possible "impurities" configurations is indicated 
by < > and He = -tt-)- The three curves for a/\ e i > 100 are 

calculations drawn from Medina et aVs work 19] and the 100 
mK data of Fig. 3 are reported by using Eq. (4). Both the 
sign and the sublinear behavior of our data are well accounted 
for by the calculations. Extrapolation from the logarithmic 
behavior at large a/A e ; gives a/A e ; ~ 5 — 10 for our sample, 
yielding an elementary dipole po ~ 10 D not that far from the 
one drawn from echo measurements. 



counts for most of the salient features of the low-T Se' r (H) 
reported in Fig. 3. 



C. Accounting for the data by using the disordered 
tunneling barrier model 

Figure 4 is drawn from Medina et aVs calculations [l9| 
and shows that < In Ao > increases with H (in Fig. 4 
H c = jx2~)j due to the additional quantum phase in- 
troduced by H. Using Eq. (4) with C = 0.03 pF and 
A e ; = 0.2 nm, the positive 5e' r (H) reported at 100 mK 
on Fig. 3 can be displayed on Fig. 4. Comparison with 
the calculations, made for large a/X e i., reveals that both 
the sign and the sublinear dependence Se' r (H) ~ H o e ± - 2 
observed on the data are compatible with disorder ef- 
fects. Moreover, using the roughly logarithmic decrease 
of < Sln(A (H)) > with a/X e i, the experimental data 
are quantitatively compatible with a ~ 5A e / ~ 1 nm, 
yielding an elementary dipole po ~ 10 D [2(]]: this is the 
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good order of magnitude with respect to the 4 D value 
drawn from pulse echo experiments Last, let us note 
that, once reported on Fig. 4, the 50 mK and 70 mK 
data lie very close to the 100 mK data, yielding the same 
8e' r {H) ~ ^jO-6±o.2 ag we vj ag the same order of magni- 
tude for a/Xei. : these low T data were not reported on 
Fig. 4 since the expression of Ctls used to derive eq. (4) 
does not hold precisely below 90 mK. 

To summarize, the microscopic model presented here 
allows to fit quantitatively the data with very rea- 
sonnable parameters. This is due to the extreme H- 
sensitivity of tunneling through a disordered barrier 
which yields non negligible effects even for H <C H c . 
Indeed, suppressing disorder and taking a = A e ; = 0.2 
nm, would lead, by standard tunneling calculations, to 
< 5\n(A (Hj) >~ -(H/H c f ~ -10~ 8 , i.e., to an un- 
detectable quadratic negative effect, at odds with experi- 
ments. 

Despite its success in accounting for our data, a limit of 
our mechanism is that it yields a monotonous H -behavior 
of 5e' r (H), i.e., it cannot account for the peaked structure 
around H = 0.03 T reported in [n). However, even if our 
scenario is based on disorder effects, it might not con- 
tradict the interaction picture used in 113. Indeed, the 
disorder explanation lies upon a Ao tma x(H) effect with a 
very large H c ~ 10 4 T while the interaction mechanism 
involves a Ao tm i n (H) effect with a very small H c ~ 10~ 2 
T (due to the ver y la rge number N of correlated TLS's). 
In the Fig. 3 of |12|, it scorns that collective tunneling 
can be viewed simply as a simultaneous tunneling of N 
elementary dipoles. If it is so, both mechanisms should 



simply add up due to their very different H c . It is well 
admitted |llj| that averaging problems in the interaction 
picture demand further progress. It should explain, e.g. 
why, in 5e' r {H) is peaked around a H c field de- 
creasing when T increases , which naively amounts to a 
counter-intuitive increase of N when T increases. 



V. CONCLUSION 

To summarize, a e r (H) dependence was shown in a 
a - Si0 2+x C y H z glass for H ~ 1 T and T < 0.4 K. 
This demanded to decrease the T uncertainty to less 
than 0.5% : beyond chemical differences between glasses, 
this might explain why previous studies reported unde- 
tectable e r (H) effects. At low T, the reported e r (H) ef- 
fects were interpreted assuming that tunneling is affected 
by some disorder within the elementary tunneling barri- 
ers of size a ~ 1 nm. This scenario might simply add up 
to the interacting one previously proposed. 
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VI. APPENDIX 
A. Proof that 77 lies in the interval [0.1; 10] 

It is argued here that, at any reported T, r\ lies in the 
interval [0.1, 10]. Defining the exponents s and t by ST(1 
T < H < 3 T) - H* and SC int .(l T < H < 3 T) - H", 
we will first assume that t — s > 0.3 and show, from 
the data of Fig. 2, that this implies 0.1 < r\ < 10 at 
any reported T. Reciprocally, assuming 0.1 < 77 < 10, 
we will show, by using the trends of Fig. 4 and Fig. 1 
that t — s > 0.3. Note that s and t are not necessarily 
integers, due to subtle disorder averaging effects , such 
as those studied in Medina et al's work (see section 
IV). 

Since we must precisely compare the relative impor- 
tances in eqs. (la)-(lb) of the first term containing 
ST(H) and the second one containing Se r (H), it is 
more convenient to translate the measured SC(H) and 
the sought Sdnt (H) into thermal equivalents by in- 
troducing ST C (H) = 5C(H)/(^) and 5T lnt . c (H) = 
SCi n t.(H)/(^). Eq. (la) thus amounts to: 

6T C (H) = 6T{H) + 5T mt . c {H), (Al) 

while defining the corresponding quantities STr(H) 
and 5Ti nt .n(H) for R measurements yields from Eq. (lb): 

5T R {H) = 5T{H) + ST mt . R (H). (A2) 

With these definitions note that the definition of a 
given in Eq. (2b) becomes a = STi nt n/5Ti nt c with still 
V = l/(a-l). 

1. Assuming t — s > 0.3 implies 0.1 < n < 10 at any 
reported T 

i) rj < 10 (i.e., a > 1.1) . Consider the magnitudes of 
the SC(H) and SR(H) bumps reported on Fig. 2 : for ex- 
ample at 50 mK, these positive bumps around H = 0.5 T 
amount to 8Tq — —6 mK and to STr ~ —7.6 mK. These 
two values are ten times larger than that derived from 
the small negative magnetoresistance of the Ru02 ther- 
mometer (see inset of Fig. 1) and, even if the inset data 
cannot be accurately trusted, these 5Tc,6Tr values are 
anyway too large to be entirely attributed to the Ru02 
thermometer |lj]. This means that both 5Ti nt .c(H) and 
STint. r(H) are negative at 50 mK, i.e., a > at 50 
mK and both 8C in t.{H) and 8R int .{H) positive in the 
Tesla range : with the above definition of s, 5Ci n t.(H) 
and SRi n t.(H) are expected to be proportionnal to the 
same H s with positive coefficients. Thus SCi nt . (H) and 



5Rint.{H) increase with H : the (quasilinear) decrease 
seen on Fig. 2 at large H thus comes from Ru02 effects, 
i.e., 6T(H) oc H l and t > s. Setting these H expansions 
in Eqs. (Al),(A2) leads to a = (Hr/HcY' 3 with H R 
the magnetic field where 8R(H) vanishes and He the cor- 
responding one for 8G(H) : at 50 mK, H R ~ 4.2 T and 
He — 2.9 T, which with the above assumption t — s > 0.3 
leads to a > 1.1, i.e., rj < 10. Finally, the same argu- 
ment can be made at T > 50 mK where H r /Hq ratios 
are larger, reaching 2.5 at 400 mK. 

ii) rj > 0.1 (i.e., a < 10). The assumption t > s + Q. 3 
implies that the low field region H < 1 T is that where 
the relative importance of 5Ti nt .c (resp. ST^.r) in the 
measured 5Tc (resp. ST R ) is the largest. In this range 
where SC(H) and 5R(H) have the same sign, the data of 
Fig. 2 correspond, for any given T, to values of STc and 
STr which differ by less than a factor 2. This means that 
STint.c and 5Ti nt .R are of the same order of magnitude 
(this might be seen also in [ll|), hence a < 10 i.e, 77 > 0.1: 
this result holds even at large H since a and 77 do not 
depend on H. 



2. Assuming 0.1 < r\ < 10 implies t — s > 0.3 

This assumption on 77 allows Fig. 3 and Fig. 4 to be 
drawn from Fig. 2, yielding s = 0.6 ± 0.2 for H > 1 
T. Besides, the data of the lower inset of Fig. 1 about 
0ruO 2 (H) effects amount to t > 1.2 for 1 T < H < 3 T 
(note that we do not rely here on the precise values of 
0RuO 2 (H) effects but only on their general trend). Hence, 
the data ensure t — s > 0.3. 



B. Ruling out spin effects at low T 

Due to the lack of microscopic information on TLS's, 
spin effects cannot be ruled out at all T. Indeed, if each 
TLS contains a single electronic spin, the non negligi- 
ble Zeeman energy SEz{H) ~ 0.7 K/T comes into play 
: each TLS becomes a four-level system. This compli- 
cates greatly the e' r calculation in the range T > T rev . 
where the relaxational (i.e., "thermodynamic") contri- 
bution plays an important role. This is why the negative 
Se' r (H) close to 400 mK has not been interpreted and the 
efforts were focused on the positive 8e' r (H) at low T. In- 
deed, well below T rev _, the TLS susceptibility is driven by 
the transitions coherently induced by £, and such quan- 
tum transitions are forbidden between states of different 
spins : TLS's can thus be separated into two independent 
subclasses of a given spin state whose susceptibilities can 
be added. Thus, for a given spin state, the asymmetry 
energy A between the two potential wells of given TLS 
is effectively independent of SEz(H). 



